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Abstract. Upon the discovery of power laws |l9][T2l[3T|, a large body of work in com- 
plex network analysis has focused on developing generative models of graphs which mim- 
ick real- world network properties such as skewed degree distributions 131], small diam- 
eter f3\ and large clustering coefficients P0ll49| . Most of these models belong either to 
the stochastic, e.g., ||9l[l4l|2l]g2), or the strategic e.g., Il6l|7iri5 30|, family of network 
formation models. 

Despite the fact that planar graphs arise in numerous real-world settings, e.g., in road 
ps| and railway maps, in printed circuits, in chemical molecules, in river networks lfT0ll43]| . 

comparably less attention has been devoted to the study of planar graph generators. In this 
^ work we analyze basic properties of Random Apollonian Networks pTl ISTI, a popular 

'""5 stochastic model which generates planar graphs with power law properties. 

(^\ Specifically, let k he a constant and Ai > A2 > .. > be the degrees of the k 

highest degree vertices. We prove that at time t, for any function / with f(t) — > +00 
as t +00, 77^ < Ai < f{t)t^/'^ and for i = 2, . . . , fc = 0(1), 7^ < A, < 



^ Ai_i — jjj^ with high probability {whp ). Then, we show that the k largest eigenvalues 

^ O ^ of the adjacency matrix of this graph satisfy At = (1 ± o(l)) AJ,^^ whp . Furthermore, we 

prove a refined upper bound on the asymptotic growth of the diameter, i.e., that whp the 
CO diameter d{Gt) at time t satisfies d{Gt) < plogt where ^ = is the unique solution 

^ greater than 1 of the equation 7^ — 1 — log 77 — log 3. Finally, we investigate other properties 

of the model. 

(N 

1. Introduction 

O In recent years, a considerable amount of research has focused on the study of graph 

^ structures arising from technological, biological and sociological systems. Graphs are 

the tool of choice in modeling such systems since the latter are typically described as a 
set of pairwise interactions. Important examples of such datasets are the Internet graph 
^ (vertices are routers, edges correspond to physical links), the Web graph (vertices are 

web pages, edges correspond to hyperlinks), social networks (vertices are humans, edges 
correspond to friendships), information networks like Facebook and Linkedin (vertices 
are accounts, edges correspond to online friendships), biological networks (vertices are 
proteins, edges correspond to protein interactions), math collaboration network (vertices 
are mathematicians, edges correspond to collaborations) and many more. 

Towards the end of the '90s, a series of papers observed that the classic models of 
random graphs introduced by Erdos and Renyi [l28l |29l and Gilbert [l35l did not explain 
the empirical properties of real- world networks ^ [TTl |3T1- Typical properties of such 
networks include skewed degree distributions [ST], large clustering coefficients [|40l |49ll 
and small average distances [|3l|49l, a phenomenon typically referred to as "small worlds". 



1991 Mathematics Subject Classification. 68R10, 68R05, 68P05. 

Key words and phrases. Complex networks; Random Apollonian Networks/Stacked Polytopes/Planar 
3-trees; Degrees; Eigenvalues; Diameter 



1 



2 



ALAN FRIEZE AND CHARALAMPOS E. TSOURAKAKIS 



Skewed degree distributions have widely been modeled as power laws: the proportion of 
vertices of a given degree follows an approximate inverse power law, i.e., the proportion 
of vertices of degree k is approximately Ck~°'. 

Understanding the properties of real world networks has attracted considerable research 
interest in the recent years [19]. A large body of work has focused on finding models 
which generate graphs which mimick real world networks. Most of the existing work 
on network formation models falls into the stochastic and the strategic categories. Kro- 
necker graphs [[42ll . the Cooper-Frieze model for the Web graph llITI . the Aiello-Chung-Lu 
model [1], Protean graphs [46] and numerous other models belong to the former cate- 
gory. The strategic approach has its origins in the work of Boorman fTSl, Aumann [[6l 
and Aumann and Myerson [jVl and a variety of models fall into this category, e.g., the 
Fabrikant-Koutsoupias-Papadimitriou model [30|. Despite the large amount of work on 
models which generate power law graphs, a considerably smaller amount of work has fo- 
cused on generative models for planar graphs. Planar graphs have been studied mainly in 
the context of transportation networks and occur in numerous real- world problems: mod- 
elling city streets [|27ll43l . crowd simulations [37J, river networks, railway and road maps 
printed circuits, chemical molecules, see also HOI and references therein. 
In this work, we prove fundamental properties of Random Apollonian Networks (RANs) 
[|52l . a popular model of planar graphs with power law properties [38J. We use the symbol 
Gt to denote the random graph at time t. The details of the model are exposed in Section[2} 
Specifically our main results are the following theorems: 

Theorem 1 (Highest k Degrees). Let Ai > A2 >■■■> be the k highest degrees of 
the Random Apollonian Network Gt at time t where k is a fixed positive integer. Also, let 
f{t) be a function such that f{t) — )■ +00 as t ^ +00. Then whp 

and for i = 2, . . . ,k 

^1/2 ^1/2 

The growing function f{t) cannot be removed, see [|32l . Using Theorem [T] we relate the 
highest k degrees and eigevalues. 

Theorem 2 (Largest k Eigenvalues). Let \i > \2 > ■ ■ ■ > be the largest k eigenvalues 
of the adjacency matrix ofGt- Then whp Aj = (1 ± o(l))-\/Ai. 

Also, we prove the following theorem for the diameter. 

Theorem 3 (Diameter). The diameter d{Gt) ofGt satisfies in probability d{Gt) < plogt 
where ^ = r] is the unique solution greater than 1 of the equation rj — 1 — log?] = log 3. 

The outline of the paper is as follows: in Section [2] we describe the model and discuss 
existing work on RANs. In Sections |3] and |4] we prove Theorems [T] and [2] respectively. 
In Section [5] we give a simple proof for the asymptotic growth of the diameter and prove 
Theorem 3] Finally in Section |6] we investigate another property of the model. Finally, in 
Section |7 we conclude by proposing several problems for future work. 
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(A) t = 1 (B) t = 2 (C) i = 3 (D) t ^ 100 



Figure 1 . Snapshots of a Random Apollonian Network (RAN) at: (a) 

t = 1 (b) t = 2 (c) t = 3 (d) t = 100. 

2. Random Apollonian Networks 

ApoUonius of Perga was a Greek geometer and astronomer noted for his writings on 
conic sections. ApoUonius introduced the problem of space filling packing of spheres 
which in its classical solution, the Apollonian packing -see [|36l for remarkable proper- 
ties of Apollonian circles-, exhibits a power law behavior. Specifically, the circle size 
distribution follows a power law with exponent of about 1.3 [16]. Apollonian Networks 
(ANs) were introduced in H and independently in [26]. Zhou et al. J52l introduced the 
Random Apollonian Networks (RANs) and gave formulae for their order, size, degree 
distribution (using heuristic arguments), clustering coefficient and diameter. High dimen- 
sional RANs were introduced in [|51l . The degree distribution of RANs was shown to 
follow a power law in [50^1521 . Other properties have also been analyzed, e.g., average 
distance of two vertices [[T2ll . properties of the connectivity profile [24]. Zhou et al. [|52|| 
proposed a simple rule that generates a random two-dimensional Apollonian networks 
with very large clustering coefficient. RANs are planar 3-trees, a special case of random 
fc-trees [39]. The general result of the degree distribution of random /c-trees was proved 
by Cooper & Uehara and Gao [[22l[33||. In RANs -in contrast to the general model of ran- 
dom k trees- the random k clique chosen at each step has never previously been selected. 
For example in the two dimensional case any chosen triangular face is being subdivided 
into three new triangular faces by connecting the incoming vertex to the vertices of the 
boundary. Random /c-trees due to their power law properties have been proposed as a 
model for complex networks, see, e.g., [|22l[34ll and references therein. Recently, a variant 
of fc-trees, namely ordered increasing k-trees has been proposed and analyzed in [[451 . 

An example of a two dimensional RAN is shown in Figure [T] The RAN generator takes 
as input a parameter Tmax equal to the number of iterations the algorithm will perform and 
runs as follows: 

• Let Go be a triangle (1, 2, 3) embedded on the plane, e.g., as an equilateral triangle. 

• fort ^Ito Tmax: 

- Sample a face Ft = (i, j, k) of the planar graph Gt-i uniformly at random. 

- Insert the vertex t + 3 inside this face (e.g., in the barycenter of the corre- 
sponding triangle) and draw the three edges (z,t + 3), (j,t + 3),(/c,t-|-3), 
e.g., as straight lines. 

It's worth pointing out why we expect this model to yield a power law degree distribu- 
tion: consider any vertex ?; at some time t. Vertex v sees a wheel graph around it (except 
for vertices 1,2,3 who see a wheel modulo one edge) as Figure |2] shows. Since we pick 
a face uniformly at random the probability that v increases its degree by 1 is (roughly) 
proportional to its degree. Therefore, we expect this process to result in a graph, which 
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Figure 2. What each vertex v sees at time t: (a) Vertices 1, 2, 3 (b) The 
other vertices {4, . . . , t + 3}. Since we pick a face uniformly at random the 
probability that the degree of v increases by one is (roughly) proportional 
to its degree. Therefore, we expect this process to generate a planar graph 
with properties similar to those generated by a preferential attachment pro- 
cess. 



is obviously planar, and which w/?p[jhas properties similar to those generated by a pref- 
erential attachment process [HI. Therefore, it should not come as a surprise that existing 
techniques developed in the context of preferential attachment models can be adapted to 
solve problems concerning the structure of Random Apollonian Networks. 

Finally, we shall make use of the following formulae for the number of vertices (n^), 
edges (mt) and faces (Ft) at time t in a RAN Gt'. 

nt = t + 3, mt = 3t + 3, Ft = 2t + 1. 

Note that a RAN is a maximal planar graph since for any planar graph nit < 3nt — 6 < 
3t + 3. 

3. Highest Degree Vertices: Proof of Theorem [I] 

We decompose the proof of Theorem [1] into several lemmas which we prove in the 
following. The proof of Theorem [T] follows from Lemmas |2} [3| |4} [5| |6] We partition 
the vertices into three sets: those added before to, before ti and after ti where = 
log log log (/(t)) and ti = log log (/(t)). We define a supernode to be a collection of 
vertices and the degree of the supernode the sum of the degrees of its vertices. 

Lemma 1. Let dt{s) denote the degree of vertex s at time t. and let a^'^^ = a{a + l) . . . (a+ 
k — 1) denote the rising factorial function. Then, for any positive integer k 

E[M.n<^-^Cj)i (1) 

Proof. We distinguish two types of vertice^ since -as it can also be seen in Figure [I]- the 
three initial vertices 1,2,3 have one face less than their degree whereas all other vertices 
have degree equal to the number of faces. 
• Case 1 g > 4: 

Note that ds{s) = 3 and that ^ > By conditioning successively we obtain 

sequence of events occurs with high probabiHty whp if lim„_j.+oo Pr [£n] — 1- 
^Despite the fact that the resuhs don't change asymptotically, we treat both cases in Lemma [l] This 
analysis will be omitted in the subsequent lemmas. 
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• Case 2s e {1,2,3}: 

Note that initially the degree of any such vertex is 2. For any A; > 

E[dtisy'^ = E[E [dtisY'^'^ldt^^is 
= E 



2t - 1 



E 
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<E[«_,W)«](i + ^)<...<(^(f)t. 



□ 



Lemma 2. The degree Xt of the supernode Vt^ of vertices added before time to is at least 
t\l^\/i whp . 

Proof. We consider a modified process y coupled with the RAN process, see also Fig- 
ure [3| Specifically, let Yt be the modified degree of the supernode in the modified process 
y which is defined as follows: for any type of insertion in the original RAN process 
-note there exist three types of insertions with respect to how the degree Xt of the su- 
pernode (black circle) gets affected, see also Figure |3]- Yt increases by 1 . We also define 
Yto- Note that Xt > Yt for all t > to. Let do = Xt„ = Ft,, = 6to + 6 and 



p* 



to 



Pr [Yt = do + r\Yt 



to 



dn 



The following technical claim is proved in the Appendix |7] 
Claim 1. 

P < 



do + r-1 
do -I 



2t + l J 



Let Ai denote the event that the supernode consisting of the first to vertices has degree 
Yt in the modified process y less than tl^'^y/i. Note that since {Xt < t}l^\^} C {Yt < 
to^^Vi} it suffices to prove that Pr Ft < tg^'^v^ = o(l). Using Claim 1 we obtain 
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' ^ - ^ V 6*0 + 5 y V 2( + 1 y 



* (6to + 5)! Ut + J ^ 
-V2t + 1/ (6to + 5)6*0+5 

„4tu+7/2+2/3t,y'' 
(6to + 5)2*''+2 



□ 



Lemma 3. No vertex added after ti has degree exceeding ^t^/^ whp . 

Proof. Let A2 denote the event that some vertex added after ti has degree exceeding 
t'^'^t^/'^. We use a union bound, a third moment argument and Lemma [l] to prove that 
Pr [^2] = 0(1). Specifically 



Pr [A,] < 5^ Pr [d,{s) > t,H'/'] = 5^ Pr [d,{sY'^ > {t,H'/'Y'^] 

S=tl S=tl 



S=tl S = tl 



□ 



Lemma 4. No vertex added before ti has degree exceeding tl^^t^^"^ whp . 

Proof. Let ^3 denote the event that some vertex added before ti has degree exceeding 
t^^t^/^. We use again a third moment argument and Lemma [l] to prove that Pr [^3] = 
oil). 
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^3/2 



s=l s=l 

tl 



3/2 

s=l s=l 



< 5!v^C(3/2)to = 0(1) 
where C(3/2) = 5-^/2 ~ 2.612. 



□ 



Lemma 5. The k highest degrees are added before ti and have degree Aj bounded by 
to ^1/2 < A, < whp. 

Proof. For the upper bound it suffices to show that Ai < t^^t^/^. This follows immedi- 
ately by Lemmas |3] and |4j The lower bound follows directly from Lemmas |2} [3] and |4} 
Assume that at most k — 1 vertices added before ti have degree exceeding the lower 
bound ^t^/^. Then the total degree of the supernode formed by the first to vertices is 
0{tl^^^/i). This contradicts Lemma 2 Finally, since each vertex s > ti has degree at 
most t'^'^y/t <^ to^t^^"^ the k highest degree vertices are added before ti whp . □ 

Lemma 6. The k highest degrees satisfy Aj < Aj_i — whp . 

Proof. Let A4 denote the event that there are two vertices among the first ti with degree 
^t^/^ and within of each other. By the definition of conditional probability and 
Lemma [3] 



Pr [A] = Pr [A|^3]Pr [^3] + Pr [Al^sjPr [^3] < Pr [^[^3] + o{l) 
it suffices to show that Pr [^41^3] = o(l). Note that by a simple union bound 



l<si<S2<ti ,__VL ^ ^ 

where pi,si,s2 = Pr [dtisi) - dt{s2) = llAs]. 

We consider two cases and we show that in both cases maxp; si,s2 = 
• Case 1 (si,s2) ^ EjGt): 

Note that at time ti there exist mt^ = 3ti + 3 < Ati edges in Gt^. 



8 



ALAN FRIEZE AND CHARALAMPOS E. TSOURAKAKIS 



to t ' 4ti 

Pz,si,.2 < X] VY[dt{si) = r ^dt{s2) = r -l\dt^{si) = di,dt^{s2) = d2] 

r=tg'^t^/2 di,d2=3 

(2) 



c(i,d2=3 
4*1 



di-l J\d2-l J\2t + l 

2tQ + ?,\idi+d2)/2 

di,d2=3 



2t + 1 



c(i,d2=3 

<t;/^tV2,2*,^2(2t;/6^1/2)8t.-2(^^)±3Nj4*l 



t \4ti 



.2t + L 

Note that we omitted the tedious calculation justifying the transition from (2) to (3) since 
calculating the upper bound of the joint probability distribution is very similar to the 
calculation of Lemma |2l 

• CASE2 (si,S2) e E{Gt) : 



Notice that in any case (si, S2) share at most two faces (which may change over time). 
Note that the two connected vertices Si, S2 share a common face only ifsi,S2 E {l,2,3}j^ 
Consider the following modified process y': whenever an incoming vertex "picks" one 
of the two common faces we don't insert it. We choose two other faces which are not 
common to si,S2 and add one vertex in each of those. Notice that the number of faces 
increases by 1 for both si, S2 as in the original process and the difference of the degrees 
remains the same. An algebraic manipulation similar to Case 1 gives the desired result. 

□ 

4. Largest Eigenvalues of the Adiacency Matrix: Proof of Theorem [2] 

In complex network analysis the spectrum of the adjacency matrix is an important as- 
pect of the network or -more generally- of the network model with several applications, 
e.g., triangle counting BSl . To the best of our knowledge the spectrum of Random Apol- 
lonian Networks has been studied only experimentally [5|. Here, we prove Theorem [2} 
We decompose the proof of the main theorem in Lemmas [7| [8| [9| \T0\ 

Having computed the highest degrees of a RAN in Section|3} eigenvalues are computed 
by adapting existing techniques [|44ll20ll , see also [[321 for a closely related analysis, and 
taking into account the special properties of the model. We decompose the proof of Theo- 
rem |2] in four lemmas. Specifically, in Lemmas |7| [8] we bound the degrees and co-degrees 
respectively. Having these bounds, we decompose the graph into a star forest and show 
in Lemmas |9]and 10 that its largest eigenvalues, which are (1 ± o{l))\^i, dominate the 
eigenvalues of the remaining graph. 



We analyze the case where si,S2 > 4. The other case is treated in the same manner. 
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We partition the vertices into three set Si, 82,83. Specifically, let 8i be the set of 
vertices added after time and at or before time tj where 

In the following we use the recursive variational characterization of eigenvalues [fTSl . 
Specifically, let Ac denote the adjacency matrix of a simple, undirected graph G and let 
Xi{G) denote the i-th largest eigenvalue of Aq. Then 

\ • AqX 

Aj(G) = mm max — - — 

s xes'.x^o X 

where 8 ranges over all (n — i + 1) dimensional subspaces of M". 

We shall use the following lemma in our proof, specifically in the proof of Lemma |9} 

Lemma 7. For any e > and any f{t) with f(t) — t- +00 ast ^ +00 the following holds 

1 _i 

whp : for all s with f(t) < s <t, for all vertices r < s, then ds{r) < s2 2. 

Proof. Set q = . We use Lemma [TJ a simple union bound and Markov's inequality to 
obtain: 



Pr 



s=f{t) 



< J] ^Pr [4(r)(^) > (si/2+V-i/2)(5)] 

s=f(t) r=l 
t s 

< ^ ^Pr > 

s=f{t) r=l 
s=f(t) r=l 

t 



2 ^ 



ft 



□ 



Lemma 8. Let 8'^ be the set of vertices in 83 which are adjacent to more than one vertex 
of 81. Then < t^e 

Proof. First, observe that when vertex s is inserted it becomes adjacent to more than one 
vertex of 81 if the face chosen by s has at least two vertices in 5*1. We call the latter 
property A and we write s & A when s satisfies it. At time ti there exist 2ti + 1 faces 
total, which consist of faces whose three vertices are all from 81. At time s >t2 there can 
be at most 6ti + 3 faces with at least two vertices in 81 since each of the original 2ti + 1 
faces can give rise to at most 3 new faces with at least two vertices in si. Consider a vertex 
s e 83, i.e., s > t2. By the above argument, Pr [\N{s) n 8i\ > 2] < f^^. Writing \8!^\ 
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as a sum of indicator variables, i.e., \S'^\ = X]l=t2 ^(^ ^ ^^'^ taking the expectation 
we obtain 



S=t2 

1 . .^ , 2t+ 1 

2t2 + 1 

By Markov's inequality: 



<(3t8+|)ln£f±4=o(^'^') 



Pr [15^1 > t^/'] < ^ = o(l). 

Therefore, we conclude that 15*31 < t^/^ w/z;? . □ 

Lemma 9. Le? F C G be the star forest consisting of edges between Si and S3 — S'^. Let 
> A2 > . . . > Afc denote the k highest degrees ofG. Then Xi{F) = (1 — o{l))y/Ai 
whp . 

Proof It suffices to show that Ai{F) = (1 - o{l))Ai{G) for i = 1, . . . , k. Note that since 
the k highest vertices are inserted before ti whp , the edges they lose are the edges between 
5*1 and the ones incident to S'^ and S2 and we know how to bound the cardinalities of all 
these sets. Specifically by Lemma [s] l^gl < t^^^ whp and by Theorem [T] the maximum 
degree in Gt,,Gt, is less than t^^^'' = t^/\ tl^^^'^ = t^/ie for = 1/16, ea = 1/32 
respectively whp . Also by Theorem jl| A j(G') > Hence, we obtain 

A,(F) > A,(G) - t'/' - t'/'' - t'/' = (1 - o(l))A,(G). 

□ 



To complete the proof of Theorem |2] it suffices to prove that Xi{H) is o{Xk{F)) where 
H = G — F. We prove this in the following lemma. The proof is based on bounding max- 
imum degree of appropriately defined subgraphs using Lemma |7] and then using standard 
inequalities from Spectral Graph Theory [18]. 

Lemma 10. Xi{H) = o(ti/^) whp . 

Proof. From Gershgorin's theorem 071 the maximum eigenvalue of any graph is bounded 
by the maximum degree. We bound the eigenvalues of H by bounding the maximum 
eigenvalues of six different induced subgraphs. Specifically, let Hi = H[Si\, Hij = 
H(Si, Sj) where H[S] is the subgraph induced by the vertex set S and H{S,T) is the 
subgraph containing only edges with one vertex is S and other in T. We use Lemma |9] to 
bound Xi{H{Si, S3)) and LemmalSlfor the other eigenvalues. We set e = 1/64. 



Ai(i7i) < Ai(i/i) <tJ/'+^ = t33/512_ 

Xi{H2) < A,{H2) < tT^%^''' = t'''/'''\ 
X,{H3)<A,{H3)<i}J'^%"' = e'"\ 
Ai(//i2)< Ai(/7i2)<tf^^ = t'''/^°'^ 

Xi{H,3)<A,{H,3)<tT^X"' = t''"'. 
Xi{Hn)<Ai{Hn)<t'/'. 
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Therefore whp we obtain 



1=1 i<j 

□ 



5. Diameter 

As shown in Appendix B of (5T\ using physicist's methodology the diameter of a RAN 
is asymptotically upper bounded by a logarithmic factor log (t). Similarly, Zhang et al. 
[ISTil show that the diameter scales logarithmically. In the following we give a simple proof 
that the diameter d{Gt) of the graph Gt created t steps is O(logt) whp and in Theorem |3] 
we give a refined upper bound for the diameter. 




Figure 4. An instance of the process for t = 2. Each face is labelled 
with its depth. 

We begin with a necessary definition for the proof of Claim [2} We define the depth of 
a face recursively. Initially we have one face a its depth is depth{a) = 1. For each new 
face j3 created by picking a face 7, we have depth{j3) = depth{^) + 1. An example is 
shown in Figure |4} where each face is labelled with its corresponding depth. 

Claim 2. The diameter d{Gt) of the graph Gt created t steps is O(logt) whp . 

Proof. Note that if k* is the maximum depth of a face then d{Gt) = 0{k*). Hence, we 
need to upper bound the depth of a given face after t rounds. Let Ft{k) be the number of 
faces of depth k at time t, then: 

l<h<t2<...<tk<tj=l ^ j=l 

By the first moment method we obtain k* = 0(log t) whp and hence d{Gt) = 0(log t) 
whp . □ 

The depth of a face can be formalized via a bijection between random ternary trees and 
RANs. Using this bijection we prove Theorem [3] which gives a refined upper bound on 
the asymptotic growth of the diameter. 

Proof. Consider the random process which starts with a single vertex tree and at every 
step picks a random leaf and adds three children to it. Let T be the resulting tree after 
t steps. There exists a natural bijection between the RAN process and this process, see 



12 



ALAN FRIEZE AND CHARALAMPOS E. TSOURAKAKIS 




Figure 5 . RANs as random ternary trees. 



The depth of T in probabihty is 



[|23l and also Figure 
solution greater than" 

BPl, pp. 284-285. Note that the diameter d{Gt 
hence the result follows. 



- log t where - 



r] is the unique 



2 -^<- - p 

of the equation r] ~ 1 — logr] = log 3, see Broutin and Devroye 

is at most twice the height of the tree and 

□ 



The above observation, i.e., the bijection between RANs and random ternary trees can- 
not be used to lower bound the diameter. A counterexample is shown in Figure |6] where 
the height of the random ternary tree can be made arbitrarily large but the diameter is 2. 
Albenque and Marckert proved in [2J that if v, u are two i.i.d. uniformly random internal 
vertices, i.e., v,u > A, then the distance d{u, v) tends to \ogn with probability 1 as the 
number of vertices n of the RAN grows to infinity. However, an exact expression of the 
asymptotic growth of the diameter to the best of our knowledge remains an open problem. 
Finally, it is worth mentioning that the diameter of the RAN grows faster asymptotically 
than the diameter of the classic preferential attachment model [,9J which whp grows as 



log* 
log log t 



see BoUobas and Riordan [lT3ll . 




Figure 6. The height of the random ternary tree cannot be used to lower 
bound the diameter. The height of the random ternary tree can be arbitrar- 
ily large but the diameter is 2. 



6. Waiting Times 

Consider the three faces created after the insertion of the first point. Let's call the 
face which gets the first vertex A and the other two faces B, C, see also Figure [T]^b). 
Also for simplicity of the notation let the configuration of Figure [TJb) correspond to time 
t = 0. Let X equal the number of steps until a new vertex picks face B or C. Clearly, 
X e {1,2,. . .}. What is the expectation E [X]? For any t > 1 
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^j^5 + 2(j-l) 2t + 3' 

Using now the fact that E [X] = Ylt^ Pr [X > t] = 1 + Pr [X > t] and substi- 
tuting we obtain that E [X] = +00. 

7. Conclusions 

In this work we studied several aspects of Random Apollonian Networks, namely the 
highest degree vertices, the spectrum, the diameter and the "waiting times". There are 
various research directions which we plan to address in future work. Indicatively, we 
report the following: (a) In contrast to the preferential attachment model, RAN have small 
separators. Specifically, by the Lipton-Tarjan theorem pT| we know the existence of a 
good separator of size 0{y/n). This is in accordance with the empirical observation that 
real world networks have good separators (TT\. We plan to investigative the minimum 
size separator in future work, (b) Consider the natural process where an adversary with 
probability a at every step deletes a vertex from the network, e.g., uniformly at random or 
in a maliscious way. Under which conditions do we obtain a connected RAN? (c) What 
is the exact asymptotic expression for the diameter of G^? 
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APPENDIX 



Proof of Claim [S 



Proof. Let r = (to = tq, . . . , r^, Xr+i = t) be a vector denoting that Yt increases by 

insertion times 

1 at Tj for i = 1, . . . , r. We upper bound the probability of this event in the follow- 
ing.Note that we consider the case where the vertices have same degree as the number 
of faces around them, Figure [2];b). The other case is analyzed in exactly the same way, 
modulo a negligible error term. 



Pr 



do + k-1 

2Tk + l 



n 

k=l 

do{do + 1) . . . {do + r - 1 

ido + r~l)\ 
{do -1)1 



r Tfc+i — 1 

n n II- 

fc=Oj=rfc + l 




k=0 j=Tk+l 



Consider now the inner sum which we upper bound using an integral: 
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< -(rfc+i + I) log(2rfc+i + l) + 
2Tk+i + l-{do + k) 



\og{2rk+i + l-(do + k))+ 
(r, + I) log {2n + 3) - + + ^^^^ + 3 - (4 + k)) 



since 



Hence we obtain ELo Eji^l+i log (l - |±f ) < ^ + ELi where 

A = (ro + i) log (2ro + 3) - ^""^ ~ ^° log (2ro + 3 - do) 

- (r.+i + i) log (2r,+i + 1) + ^'"'•^^ ^ \~ + log (2r.+i + 1 - (4 + r)) 



and 



i?. = (r. + i) log (2r, + 3) - ^^'^^ ^ K + fe) ^2r, + 3 - {do + k)) 

-{r,+ l)log{2n + l) + '^^'^^~f^^~^\ og{2n + l-{do + k-l)). 

We first upper bound the quantities for k — 1, . . . , r. By rearranging terms and 
using the identity log {1 + x) < xwe obtain 



Bk = {Tk + I) log (1 + + log (2rfc + 3) 



<i,l,o,(2.,.3)-llog(l-i|i|) 



First we rearrange terms and then we bound the term e"* by using the inequality e ^ ^^/^ > 
1 — X which is valid for < x < 1: 
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A = -(ro + i) log (1 - ^) + (r.^, + i) log (1 - ^;^^) + I log (2r^ 
-^^log (2r,+i + l-(do + r)).^ 



3v , do 



V 2i + 1 y ^ ^ V 2to + ay 

Now we upper bound the term exp [A + Yl]^=i ^k) using the above upper bounds: 



1 V " 2Tfe+3 



< 1 e 2 + 8t+4 +4(2t+l)^ _^LJ X 



2to + 3/ V 2i + 1 



i=l V 2Tfc+3 

Using the above upper bound we get that 



r ^ 

Pr < C(r, do, io, t) n [(2^-^ + 3 - (cio + + —77^) 



where 



We need to sum over all possible insertion times to bound the probability of interest p*. 
We set r'j^^Tk- for A; = 1, . . . , r. For d = o{^/t) and r = o{t^l^) we obtain: 
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P* < C(r, do, to, ^) E n [(2^^ + 3 - (rfo + k))-'/' (l + 



to+l<Tl<..<Tr<t k=l 



r 

<C{r,do,to,t) Yl l[[{24 + 3)-'/'{l + 



1 



< 



< 



< 



C{r, do,to,t) 



to-rfi+i<r(<..<r;<t-r^i 



ri 



J] (2r^ + 3)-^/^ + ^(r^ + 3/2)- 



■3/2 



to-rf- 



C{r, do,to,t) 



T ! 



d-\-r 



C{r, do,to,t) 



T ! 



(2x + 3)-V2 + _i.(^ + 3/2)-3/2 
(V2i + 3-(do + r) + 2/3y 



C(r, do,to,t) 



r do+r-3 2r 



< ~ v'-~";-"'-/ (-2^)^/2e-2 " 2t eV2«-(do+'-)+3 



r! 



< 



do + r-l\ . 2^0 + 3 . do/2 r.^ _ ^0 
do-1 7^2t + l^ 2^0 + 3^ 



to+3/2 



X 



. 2i y/2 ( do , (c?o + r)2 (4 + r)3 r(do + r - 3) , 



2r 



2 8i + 4 4(2i + l)2 



4i 



3y2t + 3- (do + r)- 



By removing the o(l) terms in the exponential and using the fact that x ^ < e we obtain 
the following bound on the probability p*. 

. ^(do + r-l\ / 2^0 + 3 yo/2 3^,„_^^^ 
^ ^^0-1 J\2t + lJ 
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